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1.  Executive  Summary 

This  report  consists  of  a  manuscript  that  describe  the  application  of  a  regional  P 
coda  wave  methodology  to  the  earthquakes  and  explosions  to  Novaya  Zernlya  and 
Nevada  test  site  events. 


1 


2 


Regional  P-coda  for  stable  estimates  of  body  wave  magnitude: 


Application  to  novaya  zemlya  and  Nevada  test  site  events 


Kevin  Mayeda 


Weston  Geophysical  Corporation,  Lexington,  MA 


Abstract 

Regional  seismic  explosion  monitoring  requires  the  discrimination  of  small  clandestine 
nuclear  explosions  from  background  earthquakes.  The  most  successful  teleseismic  discriminant, 
the  so-called  Ms:nib,  discriminant,  compares  the  long-period  surface  waves  magnitude  (Mv)  with 
the  short-period  f-based  body  wave  magnitude  (/«/>)•  There  are  many  studies  underway  to  try 
and  extend  surface  wave  magnitude  (Ms)  estimation  to  regional  distances  and  smaller 
magnitudes.  Another  problem  that  is  encountered  is  how  to  estimate  nib  so  that  the  Ms:mb 
discriminant  is  meaningful  and  consistent  with  teleseismic  measures.  For  small-to-moderate 
sized  events,  the  teleseismic  body  wave  magnitude,  mb{P),  cannot  be  effectively  measured  due  to 
low  signal-to-noise  ratio.  We  develop  a  stable  regional  alternative  based  on  the  E-coda  that 
scales  1  -to- 1  with  the  teleseismic  mb(P),  but  with  the  advantage  of  lower  variance.  Though 
mb(Lg )  and  nib(Lg- coda)  can  be  tied  to  mb(P)  for  explosions,  they  overpredict  earthquake 
magnitudes  by  -0.5-1  magnitude  units  and  degrade  the  perfonnance  of  the  Ms:mb  discriminant. 
In  contrast,  nib(P- coda)  does  not  exhibit  this  bias,  and  can  be  used  to  extend  My.mb  to  smaller 
regional  events. 

Introduction 

For  sparse  local  and  regional  seismic  networks,  a  stable  method  of  determining 
magnitude  is  necessary  for  the  development  of  discriminants,  yield  estimation,  and  detection 
threshold  curves.  Over  the  past  several  years,  the  U.S.  Department  of  Energy  (DOE)  laboratories 
have  developed  a  regional  shear-wave  coda  wave  methodology  to  obtain  the  lowest  variance 
estimate  of  the  seismic  source  spectrum  [ e.g Mayeda  et  al.,  2003;  Phillips  et  a/.,  2003;  Mayeda 
et  al.,  2007].  Unlike  traditional  magnitudes  such  as  local  magnitude  (ML)  and  teleseismic  body 
wave  magnitude  {nib,)  which  are  relative,  narrowband  measurements  that  often  have  regional 
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biases,  the  coda  methodology  provides  stable,  absolute  source  spectra  that  are  corrected  for  S-to- 
coda  transfer  function,  scattering,  inelastic  attenuation,  and  site  effects.  The  spectra  have  been 
used  to  calculate  stable  moment  estimates  (Mw),  short-period  magnitudes  (nib,  ML),  explosion 
yields,  and  radiated  seismic  energy,  ER  [Mayeda  and  Walter,  1996;  Mayeda  et  al.,  2003;  Murphy 
et  al.,  2008]  from  as  few  as  one  station.  The  coda-derived  spectra  are  calibrated  for  the 
particular  region  of  interest  and  are  in  turn  used  as  input  into  the  Magnitude  and  Distance 
Amplitude  Correction  (MDAC)  discrimination  procedure  outlined  by  Walter  and  Taylor  [2002], 

In  addition  to  MDAC’s  regional  high  frequency  discriminants,  the  traditional  teleseismic 
discriminant,  Ms'.mb,  is  currently  being  extended  to  smaller  events  at  regional  distances.  For 
example,  detailed  global  group  velocity  measurements  are  being  used  to  develop  models  for 
Rayleigh  waves  [Pasyanos  et  al.,  2003;  Stevens  et  al.,  2001;  RitzwoIIer  et  al.,  2002;  Levshin  et 
al.,  2002]  that  aid  in  the  development  of  phase-match  filters.  These  models  are  now  being 
extended  to  periods  as  short  as  7  seconds.  New  surface  wave  magnitude  formulas  [Russell, 
2006]  and  measurement  techniques  [Bonner  et  al.,  2006]  are  being  developed  that  allow 
estimates  at  these  shorter  periods  that  are  unbiased  with  respect  to  teleseismic  Ms  estimates.  The 
problem  that  we  are  experiencing  at  the  lower  magnitudes  (nib  <  4)  is  the  lack  of  unbiased  body 
wave  magnitudes  for  discrimination  purposes. 

We  could  use  Lg  and  S„  coda-derived  /«/,  estimates;  however,  this  may  actually  hinder  the 
Ms'.mb  discrimination  performance.  Though  mb  derived  from  regional  Lg  [e.g.,  Nuttli,  1973; 
Patton,  2001]  and  Lg  coda  [e.g.,  Mayeda,  1993]  have  been  calibrated  for  certain  regions,  both  are 
S-bascd  measures,  and  thus  will  be  biased  with  respect  to  earthquakes  (Figure  1).  For  example, 
the  1992  Little  Skull  Mountain  earthquake  at  the  Nevada  Test  Site  (NTS)  had  an  Mw  of  5.5,  but 
would  have  an  nib(Lg )  of  ~6.5,  whereas  the  NEIC  and  ISC  /«/,’s  for  this  event  are  5.3.  Likewise, 
if  we  calibrate  nib(Lg )  to  teleseismic  estimates  of  nib  for  earthquakes,  we  will  underestimate  the 
mb  s  for  explosions.  The  use  of  .S'-bascd  /»/,’s  in  the  traditional  Ms'.mb  discriminant  significantly 
degrades  the  discriminant’s  performance,  since  it  tends  to  move  the  explosion  and  earthquake 
populations  closer  together. 
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Figure  1.  mh(Lg- coda)  at  station  BKS  using  the  method  of  Mayeda  et  al.  [2003]  for  selected 
NTS  explosions  are  calibrated  against  the  NEIC  teleseismic  mh  (left)  and  correlate  very 
well.  However,  the  same  path  and  site  corrections  applied  to  NTS  earthquakes  results  in  a 
bias  of  ~0.5-l  magnitude  units  (right). 

Regional  m^s  have  been  calculated  based  on  the  direct  P-based  phases  such  as  P„  [ e.g 
Denny  et  al.,  1987]  and  Pg  [e.g.,  Mayeda,  unpublished  manuscript  for  the  Korean  Peninsula; 
Tibuleac  et  al.,  2001].  However,  Mayeda  [1993]  has  shown  that  these  regional  measures  have 
significant  scatter  associated  with  them,  and  thus  significant  numbers  of  recordings  would  be 
required  to  reduce  the  variance.  The  limitation  that  Bonner  et  al.  [2006]  faced  for  small  event 
analysis  using  their  Mv(VMAX)  technique  was  finding  an  unbiased  nib  magnitude.  The  objective 
of  the  current  study  is  to  find  a  more  stable  estimate  of  mb  that  will  use  regional  and  near- 
teleseismic  P-wave  data. 

Characteristics  of  Novaya  Zemlya  P-coda 

The  following  describes  preliminary  results  using  far-regional  and  teleseismic  P-coda 
waveforms  from  NORSAR,  ARCESS,  and  AWE  Blacknest  stations.  We  specifically  wanted  to 
determine  whether  P-coda  magnitudes  would  scale  with  the  teleseismic  nib  for  both  earthquakes 
and  explosions.  Second,  we  wanted  to  ascertain  if  these  P-coda  magnitudes  exhibited  less 
variance  than  their  direct  wave  counterparts.  Figure  2a  shows  array-averaged  P-coda  envelopes 
(2-3-Hz)  for  three  Novaya  Zemlya  (NZ)  explosions  (/«/~5.8)  recorded  at  NORSAR,  roughly 
2200  km  epicentral  distance,  (note:  pre-event  noise  level  differences  reflect  seasonal  variations.) 
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We  measured  relative  P-coda  envelope  amplitudes  using  the  m/,(P)  5.9  August  18,  1983 
NZ  explosion  as  a  reference  event,  though  any  event  could  have  been  used  (see  Table  1  in 
Appendix  1).  By  scaling  narrowband  envelopes  between  our  reference  event  and  the  other 
explosions  and  earthquakes,  we  were  able  to  tabulate  relative  amplitude  estimates  and  hence,  mb 
estimates.  Figure  2b  shows  P-coda-derived  /«/,  estimates  (y-axis)  relative  to  the  maximum 
likelihood  magnitude  «i/,(ML)  for  explosions  (red  squares)  and  earthquakes  (blue  triangles) 
[Lilwall  and  Marshall,  1986;  Marshall  et  ai,  1989;  Bowers,  2002].  This  regression  was  done 
using  roughly  120  seconds  of  P-coda  in  the  2-3  Hz  band  (Figure  2a).  These  preliminary  results 
are  very  promising  in  that  earthquake  nib  s  are  also  in  good  agreement  with  m/,(ML)  (Figure  2b). 
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Figure  2.  a)  Stacked  P-coda  envelopes  (2-3  Hz)  for  selected  NZ  explosions  at  NORSAR 
and  ARCESS  and  amplitudes  were  made  relative  to  the  August  18,  1983  explosion,  b)  The 
relative  mh  derived  from  the  P-coda  are  shown  for  both  explosions  and  earthquakes. 

Paths  from  NZ  to  NORSAR  are  still  at  regional  distance,  and  one  might  expect  the  P- 
wave  and  its  coda  to  be  comprised  of  waves  that  sample  the  crust  and  upper  mantle  over  a  range 
of  take-off  angles  from  the  source.  At  teleseismic  distances  however,  we  might  expect  that  the 
averaging  nature  observed  for  local  and  regional  coda  waves  to  breakdown.  At  these  distances, 
first  arriving  P- waves  are  likely  emanating  from  a  limited  range  of  take-off  angles  near  the 
bottom  of  the  focal  sphere.  To  investigate  this,  we  processed  roughly  30  NZ  explosions  recorded 
at  the  U.K.  arrays,  Eskdalmuir  in  Scotland  (EKA)  and  Yellowknife  in  Canada  (YKA)  located  at 
~30  and  44  degrees  from  NZ,  respectively. 


P-coda  vs  mb(ML) 
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Figure  3  shows  P-coda  envelopes  at  EKA  for  4  NZ  explosions  with  roughly  the  same 
magnitude  that  were  located  within  a  few  kilometers  of  each  other.  We  see  an  immediate 
discrepancy  for  the  September  24,  1979  event.  Though  it  has  the  largest  mz,(ML)  it  is  roughly  a 
factor  of  3  smaller  in  amplitude  (0.5  in  logio)  at  EKA  relative  to  the  other  three  events.  The 
direct  P-wave,  coda,  and  PcP  phase  (not  shown)  are  all  small.  In  fact,  the  EKA  station 
magnitude  for  this  event  is  also  low  relative  to  the  global  mj(ML)  estimate.  The  closest  event  is 
the  September  27,  1978  event  but  this  does  not  appear  to  be  anomalous.  Careful  inspection  of 
the  raw  data  shows  nothing  unusual  for  the  September  24th  event,  (note:  the  pre-event  noise  is 
lower  for  the  October  11,  1982  event  because  of  improvements  to  the  electronics  in  late  1979). 
We  note  that  this  event  at  NORSAR  is  in  good  agreement  with  the  m^ML )  as  well  as  at  YKA. 
Assuming  this  is  real,  then  this  suggests  a  near-source  process  such  as  focusing  directly  beneath 
this  event.  Moreover,  the  scale-length  must  be  small  since  a  nearby  event  is  not  affected.  This 
supports  the  notion  that  teleseismic  P-codas  will  not  have  the  same  averaging  properties  that 
local  and  regional  codas  exhibit. 


Figure  3.  Teleseismic  P-coda  envelopes  at  EKA  for  4  NZ  explosions  with  roughly  the 
same  magnitude  that  were  located  within  a  few  kilometers  of  each  other.  We  see  an 
immediate  discrepancy  for  the  September  24, 1979  event  (blue)  suggesting  a  break-down  in 
the  coda’s  ability  to  average  over  the  source  and  path  effects. 

Our  preliminary  findings  suggest  that  at  regional  distances  the  P-coda  can  be  used  as  a 
surrogate  for  teleseismic  mi,  for  both  earthquakes  and  explosions  based  on  the  findings  at 
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NORSAR  for  NZ  events  ( e.g .,  Figure  2b).  At  teleseismic  distances  however,  the  P-coda  appears 
to  share  the  same  radiation  pattern  as  the  direct  P- wave  and  does  not  appear  to  average  over  the 
focal  sphere  as  is  observed  for  local  and  regional  shear  waves.  Nonetheless,  the  derived  body 
wave  magnitude  mb{P- coda)  at  EKA  and  YKA  for  NZ  explosions  (not  shown)  are  in  good 
agreement  with  the  globally  averaged  results  using  direct  teleseismic  P,  though  no  improvement 
in  scatter  is  expected. 

Characteristics  of  Nevada  Test  Site  (NTS)  P-coda 

We  next  focus  on  near-regional  P-coda  from  earthquakes  and  nuclear  tests  at  the  NTS 
recorded  by  selected  regional  broadband  stations.  Using  a  single  station  at  roughly  550  km 
(BKS)  we  derived  nib(Pg- coda)  relations  for  narrow  band  envelopes  ranging  between  1  and  3  Hz. 
At  this  distance,  we  had  roughly  60  seconds  of  P-coda  before  the  direct  Lg  arrival.  As  with  the 
NZ  study,  we  made  relative  P-coda  envelope  amplitude  measurements  for  selected  earthquakes 
and  explosions  which  all  had  independent  teleseismic  estimates  of  nib  from  the  USGS  NEIC 
catalog.  Figure  4  shows  magnitude  results  from  station  BKS  operated  by  the  Berkeley 
Seismological  Laboratory.  As  found  with  NZ,  the  near-regional  P-coda  magnitudes  do  not  show 
a  bias,  in  sharp  contrast  to  mh{Lg )  and  mb{Lg- coda)  {e.g..  Figure  1). 


Figure  4.  Single  station  estimates  of  mb(P- coda)  for  both  NTS  earthquakes  and  explosions 
plotted  against  the  NEIC  teleseismic  mb. 
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Finally,  we  compare  interstation  amplitude  measurements  to  test  the  extent  to  which 
regional  P-coda  can  reduce  scatter  compared  to  the  direct  P-wavc.  Figure  5  shows  narrowband 
amplitudes  at  stations  ELK  and  KNB,  roughly  400  and  240  km  distance,  respectively. 


Figure  5.  a)  Interstation  scatter  direct  Pg  is  roughly  two  times  larger  than  P-coda  for 
the  same  NTS  events  and  regional  stations,  ELK  and  KNB. 

We  found  that  the  narrowband  regional  P-coda  amplitudes  are  roughly  two  times  smaller 
in  data  standard  deviation  than  their  direct  wave  counterparts.  In  contrast,  for  shear  wave  coda 
we  typically  observe  a  factor  of  3-to-4  improvement  [Mayeda  et  ai,  2003].  This  difference 
could  be  due  to  P-coda  being  more  forward  scattered  and  not  spatially  averaging  to  the  extent  of 
shear-wave  codas.  This  is  supported  by  array  analysis  from  regional  waveforms  by  Wagner 
[1997]. 

Conclusions 

For  small-to-moderate  sized  events,  an  unbiased,  P-based  regional  magnitude  is 
necessary  to  seamlessly  tie  to  teleseismic  estimates  of  nib  for  seismic  discrimination  and 
explosion  yield  studies.  Currently  there  is  a  debate  within  the  explosion  monitoring  community 
as  to  whether  the  explosion  and  earthquake  populations  in  the  Mynit  discriminant  merge  or  stay 
separated  at  smaller  magnitudes  (<  nib  ~3.5).  However,  due  to  limited  numbers  of  stations  for 
regional  explosion  monitoring,  direct  phase  magnitudes  such  as  m/>(P„)  and  m/,(Pg)  exhibit  high 
variance  due  to  strong  lateral  complexity  and  source  radiation  pattern.  Though  regional  shear 
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wave  magnitudes  [ e.g mi,(Lg)  and  /n/,(L,,-coda)]  can  be  tied  to  either  explosion  or  earthquake 
nib  s,  the  fact  that  these  are  shear  wave  measurements  introduces  a  significant  bias  [e.g.,  Figure 
1]  and  will  degrade  the  performance  of  the  Ms:mb  discriminant.  We  have  found  a  regional 
equivalent  to  the  teleseismic  /«/,  using  P-coda  envelopes  which  are  roughly  two  times  less 
scattered  than  their  direct  wave  counterparts  and  scales  1  -to- 1  with  teleseismic  estimates  (e.g., 
Figure  5b).  However,  at  teleseismic  distances,  we  find  evidence  that  the  averaging  properties  of 
coda  appears  to  break  down  perhaps  due  to  sampling  only  a  narrow  portion  of  the  bottom  of  the 
focal  sphere  (e.g.,  Figure  3).  Our  next  step  will  be  to  apply  the  new  P- coda  methodology  to 
other  test  sites  and  assess  the  perfonnance  of  the  Ms:mb,  discriminant  for  smaller  magnitude 
events.  In  addition  to  discrimination,  the  stable  estimation  of  explosion  yield  for  small  tamped 
events  may  benefit  from  the  use  of  the  regional  P-coda  envelope  and  studies  are  currently 
underway. 
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Appendix  1 


Table  1 
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,  1976 

02:59:57.700 
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73.398 
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15 
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-0.45 
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16 
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10 
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?? 

-1.6  - 
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(Mikhailov, 1999 

list 

of  nukes)  (Norsar  : 

report  lowered 

mb  from  4.2 

to  3.8) 
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25 
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n25 
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24 
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73.317 

54.805 

?? 

nl2 

NA 

NA 

NA 

NA 

399186 

NOV 

15 

(319) 

,  1978 

08:30:00.000 

73.4 

55.0 

3.6 

?  ? 

-1.72 

??  -1.8 

?? 

(NORSAR 

CD, 

chemical? ) 

Earthquakes : 

185081 

AUG 

01 

(213) 

,  1986 

13:56:37.800 

73.031 

56.726 

4.26 

?? 

-1.0 

-0.75 

-1.05 

-0.75 

(Marshall  et 

al.,  1989) 

399156 

JUN 

13 

(164) 
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19:22:37.900 

75.2 

56.7 

3.5 

?? 

?? 

?? 

-2 . 1 

?? (Ringdal,  1997 

399161 

JAN 

13 

(013) 

,  1996 

17:17:23.000 

75.2 

56.7 

2.4 

?? 

?? 

?? 

?? 
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(228) 
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02:10:59.910 

72 . 648 

57.352 

3.3 

?? 

-2.2 

?? 

-2.25 

?? 

(Bowers,  2002)  3.5  (Ringdal,  1998)  LLNL  envelope  screwed  up! 
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Appendix  2: 


Over  the  past  several  years,  the  DOE  labs  have  developed  a  regional  coda  wave 
methodology  to  obtain  the  lowest  variance  estimate  of  the  seismic  source  spectrum.  The  coda  is 
the  scattered  wave  train  that  arrives  after  the  direct  arrivals,  presumably  the  result  of  scattering 
from  heterogeneity  in  the  Earth.  Thus,  regional  Mw  and  mi,  estimates  derived  from  Sn  and  Lg 
coda  are  very  stable,  even  when  only  a  single  station  is  used.  However,  these  m\,  ’s  are  inherently 
biased  for  earthquakes  because  they  are  an  .S'-bascd  measurement,  and  explosions  are  relatively 
depleted  in  S- waves.  Previous  research  projects  have  used  region-specific  nib  scales  based  on 
direct  measurements  of  Pn  and  Pg  to  improve  the  Ms:nib  discrimination,  even  though  the  /«/, 
estimates  often  had  a  large  variance. 

Figure  1  shows  results  for  Nevada  Test  Site  (NTS)  explosions  recorded  at  regional 
distances.  Here  we  compare  the  inter-station  performance  between  mt{Lg),  mb(Pn)  and 
nib(Lgcoda )  from  Mayeda  (1993).  We  see  that  the  coda-based  mb  s  have  the  lowest  standard 
deviation  by  roughly  a  factor  of  4-to-5.  This  property  makes  it  ideal  for  monitoring  situations 
where  station  coverage  is  sparse. 
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Figure  1.  (from  Mayeda,  1993) 

The  next  obvious  step  to  be  implemented  in  the  coda  wave  methodology  is  the  use  of  P 
coda  for  mb  estimates.  The  following  figures  and  text  describes  the  results  of  a  two  week  visit  to 
AWE  Blacknest  where  far-regional  and  teleseismic  P-coda  were  investigated.  We  specifically 
wanted  to  know  whether  P-coda  magnitudes  would  scale  with  the  teleseismic  mi,  for  both 
earthquakes  and  explosions.  Second,  we  wanted  to  know  if  these  P-coda  magnitudes  exhibited 
less  variance  than  their  direct  wave  counterparts. 


Figure  2  below  shows  array  averaged  envelopes  (2-3-Hz)  for  two  Novaya  Zemlya  (NZ) 
explosions  (nib- 5.8)  recorded  at  NORSAR,  roughly  2200  Km  distance.  Notice  that  both  P  and  S 
codas  are  very  similar  in  character,  (note:  pre-event  noise  level  differences  reflect  seasonal 
variations.) 
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Figure  2. 

We  measured  relative  P-coda  envelope  amplitudes  using  the  October,  24,  1990  NZ 
explosion  as  a  reference  event  (Table  1).  By  scaling  narrowband  envelopes  between  our 
reference  event  and  the  other  explosions  and  earthquakes,  we  were  able  to  tabulate  relative  coda 
amplitudes.  Figure  3  below  shows  coda  envelope  amplitude  residuals  (y-axis)  relative  to  the 
maximum  likelihood  magnitude  mj(ML)  for  explosions  (red  squares)  and  earthquakes  (blue 
triangles)  (Lilwall  and  Marshall,  1986;  Marshall  et  ah,  1989;  Bowers,  2002).  This  regression 
was  done  using  roughly  100  seconds  of  P-coda  in  the  2-3-Hz  band.  These  preliminary  results 
are  very  promising  in  that  earthquake  mb  s  are  also  in  good  agreement  with  mb(ML).  This  is  in 
sharp  contrast  to  results  from  regional  nib(Lg )  and  m/,(Pgcoda)  (e.g.,  Patton,  1988;  Mayeda  1993). 
In  those  studies,  was  tied  to  explosions  at  the  Nevada  Test  Site  (NTS)  (see  Figure  1),  however 
applying  the  same  formulas  to  earthquakes  results  in  an  overestimation  of  ~1  magnitude  unit. 
For  example  the  1992  Mw 5.5  Little  Skull  mountain  earthquake  at  NTS  would  have  an  nib{Lg )  of 
~6.6. 
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Figure  3. 

Paths  from  NZ  to  NORSAR  are  still  at  regional  distance  and  one  might  expect  the  P-wave 
and  its  coda  to  be  comprised  of  waves  that  sample  the  crust  and  upper  mantle  over  a  range 
of  take-off  angles  from  the  source.  At  teleseismic  distances  however,  we  might  expect  that 
the  averaging  nature  observed  for  local  and  regional  coda  waves  to  breakdown.  At  these 
distances,  first  arriving  P-wavcs  are  likely  emanating  from  a  limited  range  of  take-off 
angles  near  the  bottom  of  the  focal  sphere.  To  investigate  this,  we  processed  roughly  30 
NZ  explosions  recorded  at  the  U.K.  arrays,  Eskdalmuir  in  Scotland  (EKA)  and  Yellowknife 
in  Canada  (YKA)  located  at  ~30  and  44  degrees  from  NZ,  respectively. 
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